Self-bound droplet of Bose and Fermi atoms in one dimension: 
Collective properties in mean-field and Tonks- Girardeau regimes 
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We investigate a dilute mixture of bosons and spin-polarized fermions in one-dimension. With 
an attractive Bose-Fermi scattering length the ground-state is a self-bound droplet, i.e. a Bose- 
Fermi bright soliton where the Bose and Fermi clouds are superimposed. We find that the quantum 
fluctuations stabilize the Bose-Fermi soliton such that the one- dimensional bright soliton exists for 
any finite attractive Bose-Fermi scattering length. We study density profile and collective excitations 
of the atomic bright soliton showing that they depend on the bosonic regime involved: mean-field 
or Tonks-Girardeau. 
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Ultracold vapors of alkali-metal atoms, like Rb, 
85 Rb, 40 K, 23 Na, 6 Li, 7 Li, etc. are now actively stud- 
ied in the regime of deep Bose and Fermi degeneracy 
0, SB 0, HI- Trapped Bose-Fermi mixtures, with Fermi 
atoms in a single hyperfine state, have been investigated 
by various authors both theoretically 
[H and experimentally QUE El 
it has been predicted that self-bound droplets, also called 
atomic bright solitons, can be formed within a mixture 
of degenerate Bose-Fermi gases provided the gases at- 
tract each other strongly enough and that there is an 
external transverse confinement 0, [2(| Hl[ . Formation 
of bright solitons in a dilute spin-polarized Fermi gas is 
prevented by Pauli repulsion. The formation of bright 
soliton in a Bose-Fermi mixture is related to the fact that 
the system can lower its energy by forming high-density 
regions (bright solitons) when the Bose-Fermi attraction 
is sufficient to overcome the Pauli repulsion among Fermi 
atoms and any possible repulsion among the Bose atoms. 
A common point of these papers [lj|, [2(| HH is that the 
Fermi cloud is three-dimensional (3D). In fact, for not 
too strong Bose-Bose repulsion the transverse width of 
the Fermi component significantly exceeds the transverse 
width of the Bose component [Is], [2(| [2l[ . 

In the strict one-dimensional (ID) regime, the Bose- 
Fermi mixture requires an appropriate theoretical de- 
scription. The exponent of the power-law which describes 
the bulk energy of a Fermi gas as a function of its density 
depends on the dimensionality (see for instance p2j). In 
addition, even at zero temperature, the ID Bose gas can 
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never be a true Bose-Einstein condensate due to phase 
fluctuations [l|, [23[. For a repulsive ID Bose gas, one 
must distinguish two regimes: a quasi-Bose-Einstein con- 
densate (BEC) regime, well described by the ID Gross- 
Pitaevskii equation with positive nonlinearity (23l |. and 
a Tonks-Girardeau (TG) regime at very low densities, 
where the ID bosons behave as ID ideal fermions [HHHI]. 
An attractive ID Bose gas is instead well described by 
the Hartree mean-field theory [2(| [27|, i.e. the ID Gross- 
Pitaevskii equation with negative nonlinearity [281 ]. 

The existence of the TG regime above has been re- 
cently experimentally confirmed .29] in a study of the ID 
degenerate 87 Rb system. In a subsequent study of this 
system [30j], the ID Bose gas in the TG regime has been 
found to possess the peculiar property of not attaining 
a thermal equilibrium even after thousands of collisions. 
This is often termed fermionization of ID bosons in the 
TG regime. It is well-known that due to Pauli princi- 
ple the spin-polarized trapped fermions do not interact 
at low temperature and hence fail to reach a thermal 
equilibrium necessary for evaporative cooling leading to 
a degenerate state. The necessary thermal equilibrium 
was attained only in Bose-Fermi [lj| [l?], SB] or Fermi- 
Fermi [31] mixtures through collision between bosons and 
fermions or between fermions in different quantum states, 
respectively. 

In this paper we consider a Bose-Fermi mixture 
strongly confined by a 2D harmonic potential in the 
transverse cylindric radial coordinate. The ensuing ef- 
fective ID system is described in the quantum hydrody- 
namical approximation, i.e. the time-dependent density 
functional approach based on real hydrodynamic vari- 
ables or complex scalar fields. Quantum hydrodynam- 
ics is very useful for the study of static and collective 
properties of a Bose-Fermi mixture and it has been used 
successfully in 3D [2(| HH, H3] for a description of bright 
and dark solitons and collapse. We investigate the ID 
mixture of bosons and spin-polarized fermions by using 
an effective ID Lagrangian 0, EI M, EH- A Gaus- 
sian variational approach is adopted to derive axial static 
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and dynamical properties of the mixture with attractive 
Bose-Fermi scattering length (a&f < 0). The solution 
of the variational scheme was found to be in satisfac- 
tory agreement with the accurate numerical solution of 
the hydrodynamic equations. We find that a self-bound 
droplet, i.e. a Bose-Fermi bright soliton, exists also for 
very small values of |a&/|. In this case the axial width of 
the Fermi component is very large while the axial width 
of the Bose component depends on the sign and magni- 
tude of the Bose-Bose scattering length a&. Remarkably, 
the TG regime is essential to preserve a localized Bose- 
Fermi soliton for very small |a&/|; in fact, for a repulsive 
Bose-Bose interaction in the quasi-BEC ID regime the 
theory predicts a minimum value of \a b f\ below which 
the mixture is uniform, i.e. fully delocalized. For large 
values of | a;,/ 1 the Bose-Fermi system is self confined in a 
very narrow region and therefore the local axial densities 
of bosons and fcrmions strongly increase. We must re- 
member, however, that above a critical axial density the 
Fermi system is no more strictly one-dimensional, and 
the same happens for repulsive bosons. 

The paper is organized as follows. In Sec. II we present 
the model used to study the degenerate Bose-Fermi sys- 
tem. Then we derive a set of coupled equations for the 
mixture starting from a Lagrangian density. In Sec. Ill, 
by using a Gaussian variational ansatz, we demonstrate 
that for an attractive Bose-Fermi interaction, the ground 
state of the model is a self-bound bright soliton (in the 
absence of an external longitudinal trap). In Sec. IV we 
present a study of the system in the quasi-BEC regime 
and in the TG regime for attractive Bose-Fermi inter- 
action. These results are further explored in Sec. V 
considering a single Fermi atom inside the Bose cloud. 
In Sec. VI we consider the problem of coupled breath- 
ing oscillations of the Bose-Fermi system and calculate 
the frequencies of these oscillations. Finally, in Sec. VII 
we present a brief summary of our investigation and dis- 
cuss the experimental conditions necessary to achieve the 
one-dimensional Bose-Fermi soliton. 



II. BOSE-FERMI LAGRANGIAN FOR ID 
HYDRODYNAMICS 

We consider a mixture of N b bosons of mass mj, and 
Nf spin-polarized fcrmions of mass m/ at zero tempera- 
ture trapped by a tight cylindrically symmetric harmonic 
potential of frequency lu± in the transverse direction. We 
assume factorization of the transverse degrees of free- 
dom. It is justified in ID confinement where, regardless 
of the longitudinal behavior or statistics, the transverse 
spatial profile is that of the single-particle ground-state 
[ToL l35l l3rl [37l HH . The transverse width of the atom dis- 
tribution is given by the characteristic harmonic length 
of the single-particle ground-state: a±j — yf%J (2rrijuj±), 
with j = b, f. The atoms have an effective ID behavior 
at zero temperature if their chemical potentials are much 
smaller than the transverse energy Kui± [ToL [35j . 



We use a hydrodynamic effective Lagrangian to study 
the static and collective properties of the ID Bose-Fermi 
mixture. In the rest of the paper all quantities are di- 
mensionless. In particular, lengths are in units of a±b, 
linear densities in units of ajj, times in units of wj 1 and 
energies in units of hujj_. The Lagrangian density C of 
the mixture reads 

C = C b + C f + C bf . (1) 
The term Cb is the bosonic Lagrangian, defined as 

C b = r b (id t + d 2 z )^ b -\i> b \ 6 G (^jt_^-V b \M 2 > (2) 

where ip b (z, t) is the hydrodynamic field of the Bose gas, 
such that n b (z,t) — \ipb(z,t)\ 2 is the ID density and 
v b (z,t) — id z ln.(ip b (z,t)/\tjj b (z,t)\) is its velocity. Here 
gb = 2ab/a±b is the scaled inter-atomic strength with ab 
the Bose-Bose scattering length. We take gb\ < 1 to 
avoid the confinement-induced resonance [3S [.Interact- 
ing bosons are one-dimensional if gj ni < 1 3^, H|| [37| • 
For x > the function G(x) is the so-called Lieb-Liniger 
function, defined as the solution of a Fredholm equation 
and such that G(x) ~ x for < x -C 1 and G(x) ~ 7r 2 /3 
for x > 1 [H. For x < we set G(x) = x [H H]. 
Vb{z) is the longitudinal external potential acting on the 
bosons. In the static case the Lagrangian density Cb 
reduces exactly to the energy functional recently intro- 
duced by Lieb, Seiringer and Yngvason (4(| . In addition, 
Cb has been successfully used to determine the collective 
oscillation of the ID Bose gas with longitudinal harmonic 
confinement [Hj]. 

The fermionic Lagrangian density Cf is given instead 

by 

C f = r f (dt + \ m dl)i> f - ^|^/| 6 - W/l 2 , (3) 

where A m = rrib/mf and ipf(z,t) is the hydrodynamic 
field of the ID spin-polarized Fermi gas, such that 
nf(z,t) — \ipf(z, t)\ 2 is the ID fermionic density and 
Vf(z,t) = iX m d z ln(ipf(z,t)/\i/jf(z,t)\) is the velocity of 
the Fermi gas. The non-interacting fermions are one- 
dimensional if (7r 2 A m /2)n 2 <C 1 [lfj. Vf(z) is the longitu- 
dinal external potential acting on fermions. In the static 
case and with Vf(z) = the Lagrangian Cf gives the 
correct energy density of a uniform and non-interacting 
ID Fermi gas. More generally, the Euler-Lagrange equa- 
tion of Cf yields the hydrodynamic equations of the ID 
Fermi gas |34| . 

Finally, the Lagrangian density Cbf of the Bose-Fermi 
interaction reads 

C b f = -g b f\^ b \ 2 \^f\ 2 , (4) 

where g b f = ^cibf/a±b is the scaled inter-atomic strength 
between bosons and fermions, with abf the Bose-Fermi 
scattering length ■ 
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Euler-Lagrange equations of the Lagrangian C provide 
the two coupled partial differential equations for ipb and 
ipf. 



id t ip b 



9b 



2*"< T "< ~ V2|^ b | 2 



V b +9bf\^f\ 1 



1>b ,(5) 



We insert the Gaussian fields ipj(z,t) into the La- 
grangian C and integrate over the spatial variable z and 
get an effective Lagrangian [42L [43j which depends on 
<Jj{t), Zj(t), fa{t), Oj{t) and their time derivatives. By 
writing the eight Euler-Lagrange equations one finds that 
the slopes fa(t) and the curvatures 0j(t) of the fields 
ij)j(z,£) can be obtained from the widths <Xj(t) and the 
center coordinates Zj (t) through the equations 



idtipf= [-A m ^ + 7r 2 A m |V/| 4 + Vf + 9bf\H 2 ] (6) 

For gbf = and < Qb < I , the first partial differential 
equation (O reduces, in the regime gb/n-b <C 1, to the fa- 
miliar mean- field ID Gross-Pitaevskii equation i.e. to 
the ID cubic nonlinear Schrodinger equation describing 
a quasi-BEC. Instead, in the regime where everywhere 
.96 /fib f , Eq. ([5l) for bosons becomes the quintic non- 
linear Schrodinger equation proposed by Kolomeisky et 
al. [33| for the dynamics of a TG gas, which is formally 
equivalent to Eq. ([6]) describing the ID noninteracting 
Fermi gas. Actually, Girardeau and Wright [4f| have 
shown that this quintic nonlinear Schrodinger equation 
overestimates the coherence in interference patterns at 
small number of particles. Nevertheless, Minguzzi et al. 
[3^ | have found that this quintic equation is quite accu- 
rate in describing the density profile and the collective 
oscillations of the fD ideal Fermi gas with longitudinal 
harmonic confinement. If we define G{x) — x for x < 
then, when g b f = and gb < 0, the Eq. (JSJ) reduces to the 
mean-field ID Gross-Pitaevskii equation with attractive 
(negative) nonlinearity, which describes quite accurately 
the attractive ID Bose gas (27l. [28j. 



III. SELF-BOUND SOLUTION: BOSE-FERMI 
BRIGHT SOLITON 

In the remaining part of the paper we set Vb (z) — 
Vf(z) = and investigate the case of a negative Bose- 
Fermi scattering length (gbf < 0). We use a time- 
dependent variational ansatz for the fields ipj(z, t) to de- 
termine the conditions under which a self-bound droplet 
of ID bosons and fermions exists. In particular we inves- 
tigate the two main regimes of ID bosons: the quasi-BEC 
regime and the TG regime. For the two fields ipj(z,t), 
with j — b, /, we use the following Gaussian ansatz 



N 



1/2 



1/4 V 2 

7T ' (J j 



exp 



2a) 



iz + iOjZ 



(J) 



where the time-dependent variational parameters are the 
longitudinal widths <Jj(t), the centers of mass Zj(t), and 
the slopes <fcj(t) and curvatures 6j(t) of the phase. It 
is obvious that the tails of the Gaussian nb(z,t) — 
\tpb(z, t)\ 2 given by Eq. |J7]) are locally in the TG regime 
but, in our terminology, a non-uniform cloud of bosons 
is in the TG regime only if everywhere its local density 
rib(z,t) satisfies the condition 9b/nb(z,t) 3> I. 



bj = -z j -2e j z j , 11, . 



(8) 



with j = b, f. The equations of motion of the parameters 
Oj(t) and zj(t) do not depend on the phase parameters 
fait) and %(t) [HSI. They are the "classical" equa- 
tions of motion of a system with effective Lagrangian 



L = T- E 



(9) 



where 



T = f (a* + 2i b 2 ) + ^ (crj + 2zj) (10) 



is the effective kinetic energy and 

E = E b + Ef + Ebf , 



(11) 



is the effective potential energy of the system. The term 
Eb involves a complicated integral of the Lieb-Linigcr 
function G(x), namely 



E = Nb * N " 



2/j2 ' ^3/2^2 

b n "b J-°o 



+00 



-°"H^"') dy - (12) 



The other two terms, Ef and Ebf, are given by 
N f X m NfTrX m 



Ef = 



2aj 3V3V 2 . ' 



and 



_ 9 bfN b Nf 

Ebf -^bT exp 



(z b - Zff 



o\ + a) 



(13) 



(14) 



We stress that the potential energy (|I1|) of the effective 
Lagrangian © can be easily obtained from ansatz (JT]) 
without including the phase parameters fait) and 0j(t). 
On the contrary, to get the kinetic energy term (jfOjl it is 
necessary to include in the ansatz the four phase param- 
eters of Eq. ([7]) [13, EH . The kinetic term is essential to 
calculate the dynamical properties of the mixture, such 
as the collective oscillations considered in Sec. VI. 

The stable stationary state of the system is found by 
minimizing the effective potential energy: 



dE 

dZj 

dE 



0, j = b,f 



J=b,f 



(15) 
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FIG. 1: (Color online). Effective potential energy E of Eq. 
(1111) as a function of soliton widths ab and 07 for parameters: 
A m = l,N b = 100, Nf = 10, g b = 0.01 and g b} = -0.2. The 
potential has a minimum at a b = 16.90 and 07 = 21.84. 
Lengths are in units of aj_ b — U / (2mbU a_) and energy in 
units of hcj±. 



Equations (|15D lead to z b = Zf and without loss of gen- 
erality we set Zj, = Zf = 0. Equations (| 16[) can then be 
rewritten as 



2tt 3 / 2 



2N[ 

n 3/2 



+ OO 



-3y q 



\2N b 



-2y 2 qi 



9b<Jb 
2N h 



dy 



-9bfN f a\ 



bf 



2N f 2 ir\ r , 



-g bf N b a\ 



(17) 
(18) 



'bf 



where a b f = ya 2 + a 2 . From Eqs. (fT7|) and (fT8| one can 

determine the widths a b and <7/ at equilibrium. Stabil- 
ity requires that the Hessian matrix of the second partial 
derivatives of the effective potential energy E(aj) is pos- 
itive definite; equivalently the Gaussian curvature 



K a = 



d 2 Ed 2 E ( d 2 E 



da 2 da 2 



\da b daf 



(19) 



must be positive. An inspection of Eqs. (jTTJ) and (fTS"]) 
shows that there are stable solutions only for g b f < 0. For 
g b f = and g b > 0, the solutions are a b — a f — +00, cor- 
responding to infinitely extended, uniform bosonic and 
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FIG. 2: (Color online). Probability density Pj(z) = 
\tpj(z)\ 2 /Nj of bosons and fermions in the self-bound mix- 
ture with g b — 0.01 and g b f = —0.2. In the upper panel: 
N b = 100 and N f = 20; in the lower panel: N b = 300 and 
Nf = 10. Solid lines: numerical results. Dashed lines: varia- 
tional results. Units as in Fig. 1. 



fermionic clouds, while for g b f — and g b < 0, the 
fermionic ID density is uniform while the bosonic cloud 
is localized with a b = 2^/2n/(\g b \N b ). For g b f < the 
Eqs. (fT7| and (fTS)) must be solved numerically. In the 
numerical calculations the Lieb-Liniger function G(x) is 
modeled by an efficient Pade approximant based on the 
exact numerical determination of G(x) [HI]. 

As the effective potential of the problem is E of Eq. 
(fTTl) , the equation (fTTJ)) together with condition implicit in 
Eq. (fTfjf minimizes the effective potential as a function of 
the two widths a b and a f . A typical plot of the potential 
for 7Y 6 = 100, N f = 10, g b = 0.01, g bf = -0.2, and 
A TO = 1 is shown in Fig. 1. Stable oscillations of the 
system are possible around the minimum. We shall study 
different features of these oscillations in the following. 

As previously stressed, bosons are one-dimensional 
under the condition g b n b -C 1, which corresponds to 
a b ^> g b N b /y/n. For a quite large width, namely for a b » 
A r b/(v / 7i : 9fc), the bosons enter in the TG regime, where 
9b /nb ^ 1- The fermions are instead one-dimensional 
under the condition (ir 2 \ rn /2)n 2 <C 1, which corresponds 

to cr/ > Nf^X m n/2. 

We shall base the present study on the Gaussian vari- 
ational approach described above, which, like any vari- 
ational approach, should be reliable for the Bose-Fermi 
ground state studied in this paper. Moreover, the an- 
alytical variational solution provides interesting phys- 
ical insight into the problem, as we shall see in the 
following. Also, the actual numerical solution of the 
full coupled dynamics is pretty complicated to imple- 
ment for all cases reported in this paper in the vari- 
ous parameter ranges. Nevertheless, we find it worth- 
while to compare the solution of the variational scheme 
with the accurate numerical solution of Eqs. ([5]) and 
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FIG. 3: (Color online). ID Self-bound Bose- Fermi droplet 
with repulsive bosons > 0). Axial widths a b and 07 of the 
bosonic and fermionic widths as a function of the Bose-Fermi 
interaction strength \gtf\, with gtj < 0. Mixture parameters: 
N B = 100, N F = 20, A m = 1 and g b = 0.01. Upper panel: 
linear-linear scale. Lower panel: log-log scale. Between the 
dot-dashed line and the dot-dot-dashed line the bosons are 
in the quasi-BEC ID regime, while above the dot-dot-dashed 
line they are in the TG regime. Fermions are in the ID regime 
above the dotted line. Units as in Fig. 1. 



(J6j) in certain cases. We solved these numerically, us- 
ing a imaginary-time integration method based on the 
finite-difference Crank-Nicholson scheme, as described in 
Ref. [Hj]. We discretize the mean- field equations using 
a time step At = 0.05 and a space step Az = 0.05, and 
z e [-L/2,L/2] with L = 2000. The boundary con- 
ditions are ipj(—L/2) = i/jj(L/2) = 0, with j = b,f. 
We start with broad Gaussians as initial wave functions. 
In the course of the imaginary-time evolution the self- 
bound mixture is quickly formed but, due to strong Pauli 
repulsion among identical spin-polarized fermions, the 
fermionic density profile extends to many hundredths of 
length's units. It is then essential to take a very large 
space interval [— L/2, L/2] of integration to see that these 
long tails of the fermionic cloud are indeed decaying to 
zero. 

In Fig. 2 we plot two sets of numerical results for 
the probability density in the quasi-BEC ID regime. 
The figure shows that, for fixed values of the interaction 
strengths g b and g b / , the axial width of the bosonic prob- 
ability density becomes larger than the fermionic one by 
reducing the number Nf of fermions and increasing the 
number N b of bosons. The figure shows that the varia- 
tional approach can be used to give a reasonable estima- 
tion of the axial widths of the two clouds. 

We now turn to a study of the Bose-Fermi system by 
using the variational approach which enables us to ex- 
plore quite easily all regimes and extract physically inter- 
esting analytical results. With the intention of illustrat- 
ing the TG regime and the quasi-BEC Gross-Pitaevskii 
regime for bosons and ID regime for fermions for a spe- 
cific set of Bose-Fermi parameters. In Fig. 3 we plot 



Gb (solid line) and u/ (dashed line) from Eqs. (fT7| and 
(fT5)l as a function of \gbf\, with gbf < for the values 
N b = 100, N f = 20, g b = 0.01, and A m = 1 of the pa- 
rameters. In the upper panel the linear-linear scale is 
employed while in the lower panel we report the same 
results on a log-log scale, to better visualize the TG 
regime and the quasi-BEC regimes. When their width 
takes values between <Jb = g b N b /y/iT ~ 0.56 (dot-dashed 
line) and a b — N b /(g b y/Tr) — 5600 ( dot-dot-dashed 
line) bosons arc in the quasi-BEC ID regime, while they 
are in the TG ID regime if Ob > 5600 (above the dot- 
dot-dashed line). Fermions are in the ID regime when 
Of > NfyJ\ m ir/2 ~ 25, i.e. above the dotted line. For 
large |<76/|, both the widths tend to zero corresponding 
to narrow soliton(s). 

IV. QUASI-BEC ID REGIME AND 
TONKS-GIRARDEAU REGIME 

The expression for the effective Lagrangian of the Bose- 
Fermi system given by Eqs. (fTTj) — ([14]) is greatly sim- 
plified in the quasi-BEC ID regime and also in the TG 
regime. In the quasi-BEC ID regime one has G(x) ~ x 
and the expression for the effective energy E of Eq. (fTTj) 
can be evaluated analytically and written as 

1 a 8 1 7 8 

where a = N b , j3 = g^ / (2^) , 7 = A m iV / [l + 
271-^/(3^3)], and 5 = \g bf \N b N f /^ with g bf = -\g bf \. 
Then at equilibrium one finds from Eq. (|16[) that cr/ and 
S can be written as functions of a b : 

"^"j^J^ ,21) 

and also 

^C +ij) ( i+ SD 3,i '■ <22) 

Eq. (|22p implies that for any finite a b and <jf one has 5 > 
3, i.e. \g bf \ > \g bf \ min = V2g b N b / {AN f ). Thus, for g b > 
the Bose-Fermi bright soliton exists only for \g b f\ > 
\gbf\min, while the system will not be bound (ab — > 00, 
07 — » 00) for \g b f\ < \gbf\min- This last situation is 
however unphysical because for a very large en, the system 
enters in the TG regime and Eq. (|2"D|) is no more valid. 

When the Bose-Bose scattering length is zero (gb = 0) 
then \gbf\min — 0. For g b < (attractive Bose-Bose 
interaction), a b is finite even for \g b f\min — 0, but it 
cannot exceed the value a b = 2^/2tt / (\g b \N b ) with the 
corresponding value of cr/ infinitely large. When \g b f\ = 
the bosonic cloud is localized (a b — 2\[2mf (\g b \N b )) while 
the fermionic cloud is fully delocalized (cr/ = +00) in the 
axial direction. 
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FIG. 4: (Color online). ID Self-bound Bose- Fermi droplet 
with attractive bosons (<#, < 0). Axial widths Ob and a/ of the 
bosonic and fermionic widths as a function of the Bose-Fermi 
interaction strength \gtf\, with gtf < 0. Mixture parameters: 
Nb = 50, Nf = 10, A m = 1 and <?& = —0.01. Fermions are in 
the ID regime above the dotted line. Units as in Fig. 1. 



Coming back to the case of a repulsive Bose-Bose scat- 
tering length (gt > 0), we observe that, after fixing gb and 
Nb, by reducing \gbf\ the bosonic width ab increases and 
the system enters in the TG regime, where G(x) ~ 7r 2 /3 
and the effective energy E of Eq. (fTTj) can also be eval- 
uated analytically as 



E = 



1 



1 (3 



1 



2 ^ 



2al 



2 a) 



(23) 



The quantities a, 7 and 8 are the same as in Eq. (f2"D|) 
while (3 = 27V h 3 7r/(3 % /3). In this case the widths ab and 
af at the equilibrium can also be determined analytically. 
They are 



a b = 



°7 = 7 



a + /3 



(a + 0) 



3/2 



3/4 



3/2 



(« + /?) 



(24) 



(25) 



and both increase as d decreases, and only when 6 = 
\gbf\ = the two widths become infinitely large. Thus 
we conclude that for any finite value of \gbf\ and gb > 
there exists a Bose-Fermi bright soliton, the existence of 
this bright soliton being guaranteed by the behavior of 
the bosonic energy term in the TG regime. 

In Fig. 4 we show the axial widths ab and af of the 
Bose-Fermi mixture with an attractive Bose-Bose scatter- 
ing length (g b < 0). We choose N b = 50, N f = 10, and 
gb = —0.01 and plot the widths as a function of |<7&/|. 
The figure shows that, as \gbf\ goes to zero, the Fermi 
width af is much larger than the Bose width ab- In fact, 



as previously shown, at \gbf \ = 0, one finds af = +00 
while ab = 2v / 27r/(|g&|^Vf)|). For large values of \gbf\, the 
Fermi width af quickly decreases and reaches the value 
af = NfyJ\ m ir/2 w 12.5 (dotted line) below which the 
Fermi system is no more strictly one-dimensional. 



V. A SINGLE FERMIONIC ATOM IN THE 
BOSE CLOUD 

The existence of the Bose-Fermi bright soliton for both 
attractive and repulsive Bose-Bose interaction, provided 
their is an attractive Bose-Fermi interaction, even van- 
ishingly small, is due to the ID effect of quantum fluctu- 
ations described by the Lieb-Liniger function G[x). To 
understand this result further we consider the case of 
a single fermionic atom (Nf = 1) interacting with the 
Bose cloud. In this case the exact stationary Schrodinger 
equation of the single-particle wave function tpf(z) of the 
fermion is given by 



(-A m 9f-| 56J ||^| 2 )^/ = e^ / 



(26) 



where e is the energy of the fermionic bound state under 
the effective potential well V e ff(z) — —\gbf\\ipb(z)\ 2 de- 
termined by the Bose cloud. We know that for gb < the 
pure ID Bose system supports a bright soliton described 

^'y /2 S echfa 

in the present case of repulsive bosons (gb > 0) attracted 
by a single fermion gbf < we adopt for the bosonic field 
ipb(z) the following ansatz 



by i>b(z) 



Guided by this result, 



^Pb(z) 



N 



1/2 



, ,-sech[ — 



(27) 



with £ the variational width of the bosonic cloud. This 
ansatz is also suggested by the existence of the analytical 
solution for the eigenvalue of the corresponding Eq. 
for the fermionic bound state. One finds 461: 



e(0 



1 

2^ 



l + \g b f\N b Z-y/l + 2\g bf \N b S 



(28) 



The energy e(£) is the sum of the kinetic energy of the 
fermionic atom and of the interaction energy between the 
fermionic atom and the bosonic cloud. 

The total energy of the system is then the effective 
energy of the Bose cloud calculated from the Lagrangian 
density ([2]) added to the energy e(£) of the Fermi atom 
and is given by 



E 



Nb_ 



Nl 



+00 



sech 6 (y)G 



ViVbSech^y) 



dy+e(0 . 
(29) 



The existence of the Bose-Fermi bright soliton depends 
on the existence of a local minimum of this total en- 
ergy E — E(£) of the system. From Eq. 
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FIG. 5: (Color online). Collective breathing frequencies wi 
and LJ2 of the self-bound Bose-Fermi droplet with attractive 
bosons (<7j, < 0). fi is the harmonic frequency of the displace- 
ment \z 1 — of the centers of mass of the Bose and Fermi 
clouds. Parameters as in Fig. 4. Frequencies are in units of 
the frequency ui± of transverse harmonic confinement. 



find that E(£) 
f's G( 9 -^- 



\ 9b ^ Nb as £ — > +oo, since for large 

9 

for any value of y. Thus, for 



•JVi,sech 2 (y)' 3 

-co, the energy goes to zero through negative val- 
ues. Again from Eq. (|29|) we see that E(£) — > +oo 
as £ — > and therefore we conclude that the function 
i?(£), being positive at the origin and vanishingly small, 
but negative, at large £'s, must possess a negative local 
minimum that is also the global minumum of the en- 
ergy. This implies that there is always a finite value of £ 
that minimizes the total energy. As previously stressed, 
this behavior strictly depends on the properties of the 
Lieb-Liniger function G(x) for large x and, as a conse- 
quence, the Bose-Fermi bright soliton exists for any neg- 
ative value gbf even for Nf = 1. 



VI. COLLECTIVE OSCILLATIONS OF THE 
BOSE-FERMI BRIGHT SOLITON 



After having established the existence of stationary 
bright solitons in a degenerate Bose-Fermi mixture, we 
study two types of small oscillations of this system 
around the stable equilibrium position and calculate their 
frequencies. The first is the stable breathing oscillation 
of the system around its mean position and the second 
describes the stable small oscillations once the centers of 
the Fermi and Bose clouds are slightly displaced with re- 
spect to each other. Given the effective Lagrangian ©, 
the problem of small oscillations is solved by expand- 
ing the kinetic energy (fl"0|) and the potential energy (fTTj) 
around the equilibrium solution up to quadratic terms. 
In this way we extract |43j the frequencies ui of the collec- 
tive breathing modes of the Bose and Fermi clouds from 



the eigenvalue equation 




N b 
N s \ m 



(30) 



where the partial derivatives must be calculated at the 
equilibrium (crf,,^), where a b and af are the Bose and 
Fermi widths obtained from Eqs. (fl~T)) and (fl~5|) . The two 
frequencies io\ and 0J2 then read 



\ 



d 2 E 



2\ m N b N s 



(31) 



where 



A 



N, 



d 2 E 
da) 



+ X m N f 



d 2 E 



AX m N b N f K G (32) 



and Kg is the Gaussian curvature of Eq. (Qj 

In addition, by taking into account the "reduced mass" 
X m N b Nf j (N b + X m Nf), one can derive the frequency O 
of harmonic oscillation of the relative distance \z b — Zf\ 
(displacement) of the two clouds from the equation 



d 2 E 



X m N b N f 

N b + x m N f " ~ d(\z b - Zf \y 



The result is 



O = 



l2\g bf \(N b + X m N f ) 



7r(er 



7 



)3/4 



(33) 



(34) 



where again a b and 07 are the Bose and Fermi widths at 
equilibrium, where z b — Zf. 

In Fig. 5 we plot the frequencies uq and L02 of the cou- 
pled breathing modes of the Bose and Fermi clouds by 
choosing the same parameters of Fig. 4, namely N b = 50, 
Nf = 10, X m = 1 and g b = —0.01. The frequencies are 
plotted as a function of the Bose-Fermi strength \g b f\. 
For \g b f\ = 0, the frequencies are decoupled: u>i is the 
frequency of the fcrmionic axial breathing mode and 0J2 
is the frequency of the bosonic axial breathing mode. 
Without a Bose-Fermi interaction the Fermi cloud is de- 
localized (07 = 00) and its breathing frequency is uq = 0, 
while the Bose cloud remains localized (due to the nega- 
tive Bose-Bose strength g b ) and its breathing frequency 
u>2 remains finite and is equal to ui2 = g b N^ / (1&tt 2 ) . Fig. 
5 shows that both the breathing frequency uq and the 
harmonic frequency Q of the displacement \zi — 22 1 start 
from zero and grow as \g b f\ increases. 



VII. CONCLUSION 

We have studied a degenerate ID Bose-Fermi mixture 
by using the quantum hydrodynamics. We find that for 
attractive Bose-Fermi interaction (g b f < 0) the ground 
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state of the system is a self-bound Bose-Fermi droplet. 
The nonexistence of a threshold in the strength of an 
attractive Bose-Fermi interaction for the formation of a 
Bose-Fermi bright soliton in one dimension is confirmed 
in the case of a single Fermi atom immersed in a de- 
generate Bose gas with repulsive Bose-Bose interaction. 
We also calculate the frequencies of stable oscillation of 
the Bose-Fermi bright soliton. Such a Bose-Fermi bright 
soliton is similar to a recently studied Bose-Bose bright 
soliton bound through an attractive interspecies interac- 
tion H3. 

In view of the recent experimental studies of a degener- 
ate ID 87 Rb gas [29|,[3(| and the successful identification 
of the quasi-BEC and TG regime in it and the observa- 
tion of degenerate Bose-Fermi mixture in 6 Li- 7 Li HS[l2|, 
40 K- 87 Rb [13, 6 Li- 23 Na 18]] etc. by different groups, the 
experimental realization of a Bose-Fermi bright soliton 
seems possible with present technology. The most attrac- 
tive procedure seems to make use of an experimentally 
observed Feshbach resonance |48| in a Bose-Fermi mix- 



ture. The ID Bose-Fermi mixture must be created in an 
axial harmonic trap and the Bose-Fermi interaction must 
be turned from repulsive to attractive by manipulating a 
background magnetic field. At the same time the axial 
harmonic trap on the system should be removed. Upon 
removal of the axial trap, the result is the formation of 
a single or a train of bright solitons as in the experiment 
with the degenerate Bose system of 7 Li atoms [17| or 
as in a numerical simulation in a degenerate Bose-Fermi 
mixture [T^ |. By choosing numbers of atoms and inter- 
atomic strengths as suggested in the present paper, one 
obtains a single Bose-Fermi bright soliton and can study 
its static and dynamical properties. 
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